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A fully relativistic description of the electron energy loss and the induced radiation emission in 
arbitrary arrays of non-overlapping metallic cylinders is presented in terms of the multiple scattering 
method on the basis of vector cylindrical waves. Numerical analysis is given for dilute and dense 
arrays of Aluminum cylinders with a nanoscale diameter. The results of the electron energy loss 
spectrum are well correlated with the dispersion relation of coupled surface plasmon polaritons, and 
can be interpreted with an effective medium approximation when the electron runs inside the arrays. 
In addition, the cavity modes localized in the grooves between the cylinders can affect strongly the 
electron energy loss spectrum. 
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I. INTRODUCTION 

Recently, much interest has been attracted in opti- 
cal properties of composite materials such as cluster of 
metallic nano-particlesi, photonic crystal 2 ^, and left- 
handed material^*&. They have the potential ability to 
enhance various optical processes with their rich spec- 
trum, and will be a key component of future optoelec- 
tronics devices. So far, their properties have been investi- 
gated mainly with far-field optical equipment. However, 
near-field responses of them are of great importance, be- 
cause much information can be attained from the evanes- 
cent light involved in composite materials. 

The Electron energy loss (EEL) spectroscopy in a scan- 
ning transmission electron microscope is a useful tool to 
investigate both surface and bulk excitations of samples. 
In a typical EEL experiment the kinetic energy of elec- 
tron is on the order of 100 [keV], and a low-energy part 
of the EEL spectrum is related to a collective excitation 
such as surface plasmon. Thus, the spectrum is well inter- 
preted in terms of a classical macroscopic theory on the 
basis of an effective dielectric function, taking account of 
geometry and anisotropy of samples. For instance, the 
EEL of a multi-wall fullercnc can be well explained with 
an effective medium in which the dielectric constant of 
the normal direction to the multi-walls is different from 
that of the other direction^ 

There are extensive references on the EEL spectrum for 
various geometrical object s&2*i°*ii, though most works 
neglect the retardation effect. As for the composite ma- 
terials, Maxwell-Garnett(MG)-type approximations have 
been widely used. In the last decade a remarkable 
progress has been made on this subject. Pendry et 
al. proposed a new theoretical method for the EEL 
of periodically arranged nano-structuresi^ in terms of 
the real-space transfer matrix of the electro-magnetic 
wave^ Garcia de Abajo et al. developed a multiple- 
scattering method for clusters of nano-particles, on the 
basis of vector spherical wavesi^i^ and of the boundary 
element method4S*ii These methods are very powerful 
and in principle can be adapted to various complex ge- 
ometries. However, owing to the spatial discretization 



used in the transfer matrix method and in the boundary 
element method, the accuracy of the results reduces to 
some extent at high frequencies. 

Here, we present a fully relativistic multiple-scattering 
approach focused on arrays of non-overlapping cylinders 
by using vector cylindrical waves. Though at present the 
method works only for the cylinders with isotropic di- 
electric functions, in principle it can work for anisotropic 
cylinders such as multi-wall carbon nano-tubesi^. With 
this method the light scattering, as well as the EEL 
and the induced radiation emission in arbitrary arrays 
of cylinders can be treated in a unified manner by solv- 
ing exactly the multiple-scattering equation. In addition, 
since the method utilizes the Fourier decomposition along 
the cylindrical axis, the local optical responses of a cer- 
tain wave vector along the axis can be obtained. This 
property is feasible in order to see the effects of the lo- 
calized electro-magnetic modes of cylinders in detail. 

A certain amount of the EEL in composite materi- 
als is caused by the induced radiation emission. As 
for a metallic nano-particle, the phenomena is known as 
the surface plasmon radiation emission^ By arranging 
nano-particles periodically, so-called Smith-Purcell(SP) 
radiation22i2i takes place. One of the authors (K.O.) 
and his collaborators have studied the SP radiation in the 
photonic crystals composed of dielectric spheres in detail, 
where a notable enhancement of the intensity of the SP 
radiation takes place owing to the singular state density 
of photon in the photonic crystals ,22^24 However, to the 
best of our knowledge, none has been reported concern- 
ing the radiation emission from arrays of cylinders. Thus, 
the quantitative evaluation of the radiation emission in 
the arrays is another theme of the paper. 

On the other hand, the determination of the effective 
dielectric function of composite materials is still valuable 
at low frequencies. It gives a concise explanation of the 
absorption spectrum as well as the EEL spectrum of the 
composite materials in bulk, 2 - Moreover, the determina- 
tion is an important issue in the field of left-handed ma- 
terial, which has negative permittivity and permeability 
simultaneously. Several methods to determine them have 
been proposed by many authors ,2§i2L 2 £ Here, we propose 
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an alternative method by using the scattering matrix of 
semi-infinite photonic crystals and compare the EEL in 
the effective medium with that of photonic crystal. 

This paper deals mainly with the formalism of our 
method as well as the numerical analysis on clusters 
of Aluminum cylinders with a diameter of a few nano- 
meters, bearing carbon nano-tube arrays in mind. In the 
paper II we will deal with a metallic photonic crystal 
whose lattice constant is comparable with or larger than 
the plasma wavelength of the constituent metallic cylin- 
ders. Since it is very bulky to present a comprehensive 
analysis on the EEL and the induced radiation emission 
in a single paper, we should discuss the above topics sep- 
arately. In a nano-structure, which is discussed in this 
paper, the relevant range in wavelength, which is near the 
surface plasma wavelength of metal, is much larger than 
the diameter and than the pitch of the structure. As a 
result, effects of usual photonic bands, which come from 
zone foldings and degeneracy shifts on Bragg planes, do 
not appear in the EEL spectrum and the induced radi- 
ation emission spectrum at relevant frequencies. In par- 
ticular there is no remarkable feature in the SP radiation 
emission spectrum. Instead, unusual photonic bands of 
coupled surface plasmon polaritons have a strong influ- 
ence on the EEL spectrum at very low frequencies. More- 
over, in the relevant frequency range an effective medium 
approximation can be reasonably adapted to the struc- 
ture. As for the structure which will be discussed in the 
paper II, we will show that effects of both the usual pho- 
tonic bands are very pronounced for the EEL and the 
induced radiation emission spectra. 

The paper is organized as follows. In Sec. II we briefly 
summarize the vector cylindrical wave formalism. Using 
the formalism the dispersion relation of the surface plas- 
mon polariton in an isolated metallic cylinder is obtained. 
The multiple scattering method is adapted to the EEL 
and the induced radiation emission in clusters of metal- 
lic cylinders in Sec. III. In Sec. IV the expression of the 
EEL in metallic photonic crystals is derived. The nu- 
merical results of the EEL spectrum are compared both 
with those of the isolated cylinder and with those of an 
effective homogeneous medium. Finally, we summarize 
the results. 



II. SURFACE PLASMON POLARITON IN AN 
ISOLATED CYLINDER 

An infinitely long metallic cylinder with a circular cross 
section can support an electro-magnetic surface-localized 
mode that is called surface plasmon polariton(SPP). The 
SPP is characterized by an angular momentum I, a wave 
number k z , and an angular frequency to, owing to the ro- 



tation invariance with respect to the cylindrical axis, the 
translational invariance along the axis, and the transla- 
tional invariance of time, respectively. Before considering 
its dispersion relation, it is valuable to note some formu- 
las of the light scattering by an isolated cylinder using 
vector cylindrical waves>22, From now on, we take the 
cylindrical axis to the z-axis. Assume that a monotonic 
incident wave with a momentum k z along the cylindrical 
axis is scattered by an isolated cylinder with a dielectric 
function e a (uj) and a radius r embedded in a host with a 
permittivity Eb- Throughout the paper Eb is taken to be 
1 in numerical calculations, though we keep it unspeci- 
fied in equations. The incident wave can be written as a 
superposition of vector cylindrical waves: 
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where c is the speed of light in vacuum, (p, 9, z) is the 
cylindrical coordinate, z is the unit vector along the z- 
direction, and V|| is the differential operator with respect 
to (p,0). Throughout the paper, we take the follow- 
ing convention of the square root of a complex number: 
Im^w > for Im(u>) > 0. Sometimes, we call the M(N)- 
field the P(S) polarization. At k z = the M(N)-field 
corresponds to the TE(TM) polarization. 

A metallic or dielectric cylinder scatters light irrespec- 
tive of whether the light is evanescent or not. By impos- 
ing the boundary condition of Maxwell's equation, we 
can solve the scattering problem exactly. The induced 
wave scattered by the cylinder is given by 
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Here, Hi is the Hankel function of first kind and 

is the t-matrix of the cylinder, its analytical expression 

being 
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It should be emphasized that except for k z = or / = 0, 
the M(N) polarization mixes with the N(M) polarization 
in the induced wave. This property is quite distinct from 
the case of an isolated sphere, where the polarization mix- 
ing does not take place. In the light scattering by cylin- 
der there is another distinct point from that by sphere. 
Since Af, is the wave number in the (x, y) plane, the in- 
duced wave is evanescent when Af, is pure imaginary. On 
the other hand the induced wave by the sphere behaves 
as /i;((#,|x|), where hi is the spherical Hankel function of 
first kind. This means that induced wave by the sphere 
is always real and has a net flux at |x| = oo irrespective 
of whether the incident wave is evanescent or not. 

Next, we focus on the SPP mode in an isolated metallic 
cylinder. The SPP mode is a real eigenmode in the cylin- 
der and can exist without the incident light. Therefore, 
taking account of Eq.© the equation that determines 
the dispersion relation of the SPP is given by 
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Like as the SPP on a flat metal/air interface, whose dis- 
persion relation is given by the pole of the interface S- 
matrix between the metal and air, the dispersion relation 
of the SPP in a metallic cylinder is given by the pole of 
the t-matrix. The above equation leads to the following 
secular equation: 
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for the SPP mode with angular momentum /. Strictly 
speaking, we must say that the SPP mode is referred 
to the mode with pure imaginary A a , when the cylinder 
is of metal. The secular equation also determines the 
dispersion relation of the guided modes, in which A Q is 
real and positive, when the cylinder is of dielectric. 

For simplicity, from now on, we restrict our consider- 
ation to an Aluminum cylinder with diameter 2.5[nm], 
though the formalism presented in the paper can be 
adapted to any types of cylinders as long as an isotropic 
dielectric function of the cylinder is concerned. The di- 
electric function of the Aluminum can be approximated 
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FIG. 1: The dispersion relation of the surface plasmon polari- 
ton of the Aluminum cylinder with diameter 2.5[nm] in air. 
The loss-less Drude dielectric function with Tiuj p = 15[eV] was 
assumed for the Aluminum cylinder. The bold line is the light 
line. 



well with the Drude formula: 
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uj{uj + irj) ' 



(13) 



Tiujp and hr\ being 15 [eV] and l[eV], respectively. We fo- 
cus on this cylinder embedded in air [eh = 1) in numeri- 
cal calculations of the paper. The small diameter of the 
cylinder is comparable with that of a carbon nano-tube 
and the plasma frequency is 3.67 x 10 15 [Hz] and 82.7[nm] 
in terms of wavelength, which corresponds to ultra-violet 
light. Though the cylinder has a nano-scale diameter, it 
seems still reasonable to apply the macroscopic dielectric 
function to the cylinder. 

In Fig.l the dispersion relation of the SPP modes of 
the Aluminum cylinder is shown. In this case the SPP 
mode with I — 0, which is the solution of dQ NN = 0, 
is quite distinct from those with I > 1. The mode with 
I = ends at u) — and has positive slope in the (ui, k z ) 
space, whereas the other modes end at finite lo on the 
light line and have negative slopes. As I increases, the 
dispersion curves approach to lj = ui p /y/2. Though the 
dispersion curves end on the light line, they can be ex- 
trapolated inside the light-cone. With this extrapolation 
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we can understand how the SPP resonance occurs when 
real light is incident on the cylinder. The resonance also 
affects significantly the induced radiation emission when 
a charged particle passes near the cylinder. 



III. ELECTRON ENERGY LOSS IN CLUSTER 
OF CYLINDERS 

The formalism presented above serves to describe the 
EEL and the induced radiation emission when a charged 
particle runs near a cluster of cylinders. In this method 
the retardation effect is fully taken into account, though 
the recoil of the charged particle is neglected. In a typ- 
ical EEL experiment the kinetic energy of the electron 
is about 100 [keV], whereas the total energy loss is less 
than 0.01% of it, so that the assumption of neglecting 
the recoil is fairly justified. 

Let us consider a charged point particle with a charge 
e and a velocity v runs near the cluster composed of non- 
overlapping metallic cylinders aligned in the z-direction, 
whose positions are specified by two-dimensional (x and 
y) vector x a (a — 1,2, ...,N), N being the number of the 
cylinders. For simplicity, we assume that the charged 
particle does not penetrate any cylinders in the clus- 
ter and the trajectory of the particle is perpendicular 
to the cylindrical axes. Therefore, without losing gen- 
erality, the position of the particle at time t is taken to 
x t = (vt, yo, zq) in the Cartesian coordinate. As is known 
well, a running charged particle is accompanied by the 
electro-magnetic wave that is given by 
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in the time-Fourier component. Here, the superscript of 
e is referred to the sign of y — y$. If the light velocity 
in the background medium is smaller than v, the electro- 
magnetic wave becomes real, yielding the Cerenkov radia- 
tion. Here, we restrict ourselves to the region v < c/y/Fb, 
so that the electro-magnetic wave is evanescent (7 is pure 
imaginary) . 

The above expression can be transformed into a linear 
combination of the vector cylindrical waves centered at 
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Here, 0k± is the argument of two-dimensional vector 
= (uj/v, ±7) for a real 7. In the case of a pure 
imaginary 7 we must define 9k± as 
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Moreover, the ± in Ea. (|19|l and (|20|l corresponds to the 

sign of yo — Va an( i ipfa ^ s ^ ne multi-pole coefficient of the 
incident wave for cylinder a. We should note that the k z 
integration is involved in Ea. (|17fl . However, since k z is a 
conserved quantity for the cluster, each cylindrical wave 
with fixed k z is independently scattered by the cluster. 

The incident evanescent wave is multiply scattered 
in the cluster of the cylinders. By using the multiple- 
scattering method, the induced radiation field is self- 
consistently determined as follows^ If we focus on cylin- 
der a, the incident wave consists of the direct term 
fEa.Q17[l) plus the sum of the induced wave scattered 
by another cylinder a'. Therefore, the induced wave 
from cylinder a is obtained by multiplying the t-matrix 
of cylinder a to the multi-pole components of the inci- 
dent wave re-expanded around x = x a . As a result, the 
self-consistent induced wave is determined as 
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Here, Gi a ^> a i is the propagator from cylinder a' to a, 
and in its expression p aa i and 6 aa > are the magnitude 
and the argument of x a — x Q / , respectively. 

As was mentioned in Sec. II, the induced wave con- 
tains propagating components with real Xb, which have 
a net flux at p — * 00. This implies that a radiation emis- 
sion takes place when a charged particle passes near the 
cluster. As a result, the charged particle losses its en- 
ergy via the emission. In addition, if the cylinders are 
lossy having positive imaginary part in E a (oS), a part of 
the energy is absorbed in the cylinders. The total en- 
ergy loss is then the sum of the radiation emission and 
the absorption. Qualitatively, the loss can be calculated 
with the exerted force by the induced field reacting on 
the particle. Thus, its expression is given by 
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per unit angular frequency. The integral over t yields 
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Here, ± in the above equation is referred as the sign of 
2/0 - 2/a- 

On the other hand the net flux of the induced radiation 
emission is obtained by 

PemM = Hm \ { dzdOp (29) 

xRe [(E ind (x, to))* x H ind (x, u)] ■ p. (30) 

per unit angular frequency. Using the asymptotic form 
of the induced field ( - Ea. (|22(l ). the net flux turns out to 
be 
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The cutoff of the k z integral comes from the fact that 
the scattered wave becomes evanescent when \k z \ > qb- 
If there is no absorption in the cylinder, the energy loss 
must be equal to the net flux of the radiation emission, 
that is, 



P e i(uj,k z ) = P em (uj,k z ). 



(34) 



This equality serves as a criterion of the correctness and 
the convergence of numerically calculated EEL spectrum. 

In the case of an isolated cylinder, the integration over 
9 in Ea. (|32|l can be performed analytically, yielding 
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First of all, let's consider the EEL in the isolated Alu- 
minum cylinder. The EEL spectrum of the cylinder has 
a sequence of peaks at the frequencies of the SPP modes. 
This is because the t-matrix has a SPP pole on the real 
axis in the complex plane of u> as can be seen in Eq. J7J 
and (|12|l . To see the correspondence of the dispersion re- 
lation of the SPP, Fig. 2 shows the integrand P e z(w, k z ) of 
the EEL spectrum, changing the k z value from to 10 in 
units of LOp/c. Here, the velocity and the impact param- 
eter of the charged particle were taken to 0.4c and 2r, 
respectively. In addition to the main EEL peak around 
uj = {0 v j\p2, another EEL peak appears much below 
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FIG. 2: The electron energy loss spectrum in an isolated Alu- 
minum cylinder with diameter 2.5[nm] for several k z values 
are shown. The inverse relaxation time rj of the Drude di- 
electric function was taken to be hrj — l[eV]. The impact 
parameter, i.e., the distance between the cylindrical axis and 
t 
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FIG. 3: The integrated EEL and the radiation emission spec- 
tra in the isolated Aluminum cylinder. The same parameters 
as for Fig|2|were used. 



/V2. This peak comes from the SPP mode with 
I = 0, as can be understood clearly by comparing with 
the dispersion relation of the SPP mode with 1 = 0, As 
a general tendency, P e i(uj,k z ) decreases with increasing 
\k z \- 

Next, we consider how the imaginary part affects the 
percentage of the radiation emission in the EEL. To 
this end, we show in Fig[3] the integrated EEL spec- 
trum P e i(uj) and the integrated radiation emission spec- 
trum P em (uj) at 77=0, l[eV]. At 77 = these must coin- 
cide. As can be seen, when the imaginary part in e a 
is introduced, the radiation emission spectrum is almost 
unchanged from that in the loss-less cylinder except for 
the SPP frequency region. However, the EEL spectrum 
receives a drastic change. According to the figure, the 
EEL is dominated by the absorption in the cylinder un- 
der study. Therefore, the efficiency for converting the 
kinetic energy of the charged particle to the radiation 
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FIG. 4: The EEL spectra in the two identical Aluminum 
cylinders with various spatial arrangements shown in the in- 
sets. k z = was assumed. See text for the other parameters. 



via surface plasmon is very low. 

Next, we consider how the spectra change when an- 
other cylinder is added. In this case, the multiple- 
scattering of the induced radiation significantly affects 
the spectra. As we will see, almost touched metal- 
lic cylinders causes a drastic change in the spectra as 
well as in the near-field configuration. This phenom- 
ena is closely related to the surface-enhanced Raman 
scattering, 30 - 31 32 ' 33 in which the intensity of the induced 
electro-magnetic field is enhanced more than thousand 
times as large as that of the incident intensity. As an ex- 
ample, we explored the EEL for various spatial arrange- 
ments of the two identical Aluminum cylinders. Fig0| 
shows the EEL spectra P e i(uj,k z ) with k z — corre- 
sponding to the arrangements shown in the insets. In (a) 
and (b) the two-cylinders are well-separated, the distance 
between two cylindrical axes being 4r (r = 1.25[nm]), 
whereas in (c) and (d) they are very close to each other, 
the distance being 2.16r. The distance between the two 
cylindrical axes and the trajectory of the charged parti- 
cle is 2r in (a) and (b), and is I.08r in (c) and (d). As 
expected, if the two cylinders are well separated along 
the trajectory (case (a)), the EEL spectrum per cylinder 
has a single peak near u> = lu p / \/2, though an asymme- 
try of the peak is observed. In case (b) the two cylinders 
are separated along the normal direction of the trajec- 
tory. This geometry still allows a single peak in the 
EEL spectrum at ui ~ ui p / y/2. Since there is the par- 
ity symmetry with respect to the trajectory, No electro- 
magnetic modes with the odd parity is involved in case 
(b). On the other hand, if the two cylinders are very 
close to each other, several loss peaks appear at the fre- 
quencies far from ui p /y/2. In particular, the spectrum 
of the case (d) has two marked peaks at lo/lo p ~ 0.52 
and 0.84. The latter peak is shared also by case (c). 
The peak at uj/lu p ~ 0.52 is related to the weakly res- 
onant cavity modes localized in the groove between the 
two cylinders. To confirm this, we show the field inten- 
sity |if z (x)| 2 at the peak frequency in FigJSJIn the figure 
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FIG. 5: The field intensity \H z (x)\ 2 induced by a running 
charged particle whose velocity is 0.4c is plotted at (k z ,uj) = 
(0, 0.523tJp), where a peak is observed in the EEL spectrum. 
The solid line stands for the edges of the two cylinders. The 
trajectory of the charged particle is at y = — 1.35[nm]. 

we can see clearly that the peak is caused by the cavity 
mode in the groove. As for the peak at lo/uo p ~ 0.84 
we couldn't find a clear evidence of the resonant cavity 
mode, as the field intensity at the peak frequency has an- 
other local maximum at a boundary of the cylinder beside 
that in the grooves. This feature suggests that this peak 
to be caused by a strong mixing of a cavity mode and 
the SPP modes. If the electron runs across the groove 
of the almost touched cylinders (case (c)), the peak at 
lj ~ 0.52w p disappears owing to the symmetry mismatch 
and the peak at u> ~ 0.84w p receives an enhancement. 

Regarding the radiation emission spectrum in the two 
cylinders, its features are more or less similar to those 
in FigQ] though the magnitude of P em (uu,k z ) is much 
smaller than that of P e i(u), k z ) at k z = 0. 



IV. ELECTRON RUNNING OUTSIDE 
PHOTONIC CRYSTAL 

A photonic crystal that consists of a periodic array 
of metallic cylinders has a rich spectrum in it, includ- 
ing infinite SPP bands for the TE polarizatio n 34 ! 35 and a 
low-frequency cutoff for the TM polarizationi^SiSi Com- 
bining these properties with a running charged particle, 
the system can react as a novel light emitter. In fact, 
when a charged particle passes near the photonic crystal, 
it induces the emission of real photon as was first pointed 
out by Smith and Purcell for a metallic grating, 20 In the 
photonic crystal, this phenomena can be interpreted in 
two ways. One interpretation is as follows. The inci- 
dent evanescent wave from the charged particle acquires 
an Umklapp momentum transfer in the photonic crys- 
tal, thereby coming into the light cone in the (u>,k\\) 
space, fey being the wave vector parallel to the bound- 
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FIG. 6: The phase diagram of the induced radiation emis- 
sion in a periodic array of cylinders with lattice constant a 
is shown. The Smith-Purcell radiation is possible only in the 
shaded region of the (k z ,ui) plane. The velocity of the charged 
particle was taken to be v — 0.4c. The solid line is the light 
line. 



ary of the photonic crystal. As a result, real photon is 
emitted from the photonic crystal. The other interpreta- 
tion is to regard the phenomena as a coherent radiation 
emission from different cylinders. Though these two in- 
terpretations are equivalent, the two points of view give 
us a deep insight of the Smith-Purcell(SP) radiation in 
the photonic crystal. 

In the case of an isolated cylinder, the induced radi- 
ation emission is possible when q 2 > k 2 , i.e., when the 
emitted light is inside the light cone of the (w, k z ) space. 
On the other hand, the condition that the evanescent 
light turns out to be a real photon via the Umklapp mo- 
mentum transfer in the photonic crystal is given by 



q 2 b - k 2 z - (- - n— ) 2 >0, 
v a 



(36) 



a being the pitch of the photonic crystal. Therefore, the 
allowed frequency range of the SP radiation is 



,~ 2 n*e b 
u- < uj < uu +1 k z < — 



(37) 
(38) 



where Co = ujcl/2-kc and k z — k z a/2Tt. If the above con- 
dition is not satisfied in the light cone, a destructive in- 
terference among the induced radiation fields from dif- 
ferent cylinders occurs, leading to the prohibition of the 
radiation emission. This also implies that if there is no 
absorption in the photonic crystal, the concerned range 
in the (ui, k z ) space does not contribute to the EEL in 
the photonic crystal. The phase diagram of the radia- 
tion emission at v = 0.4c is shown in Fig|H| In Fig^lthe 
SP radiation is allowed in the shaded regions. The al- 
lowed region becomes narrow in the (oj, k z ) space as the 
velocity of the particle decreases. Beside, in the allowed 



region of the SP radiation the photonic density of state 
is in general singular including the Van Hove singularity. 
Therefore, we may expect a quite rich spectrum of the 
SP radiation in the photonic crystals. 

The rich spectrum is not limited in the SP radiation. 
We should mention that the photonic band structure ex- 
ists also outside the shaded region of Fig|SJ As long as 
the imaginary part in e r is non-zero, the absorption of 
the induced radiation is inevitable in the photonic crys- 
tal. This causes the EEL outside the shaded region, and 
the EEL is affected significantly by the photonic band 
structure therein. Inside the shaded region the EEL con- 
sists of the radiation emission and the absorption, and 
these are independent physical observables. 

Like as the case of cluster of cylinders, the EEL and the 
SP radiation in the photonic crystals can be treated in a 
unified framework with the multiple-scattering method 
on the basis of vector cylindrical waves. Pendry and 
Martin-Moreno presented for the first time an unified 
framework in terms of the transfer matrix method to ar- 
gue the EEL and the SP radiation in photonic crystals^, 
The method is based on a discretization of Maxwell's 
equation on a spatial mesh in order to obtain the real- 
space transfer matrix. 1 ?. They found that the scattering 
matrix, which is obtained from the transfer matrix, is 
directly related to the EEL and SP radiation spectra. 
Their algorithm is easily adapted to the two-dimensional 
counterparliSLiSiiS of the layer Korringa-Kohn-Rostoker- 
Ohtaka (KKRO) method 40 i 41 i 42 i 43 i 44 i 45 , which is a gen- 
eralization of the multiple scattering method to peri- 
odic systems. The layer-KKRO method has very high 
accuracy for the photonic crystal under consideration. 
We should note that the three-dimensional layer-KKRO 
method was already used for discussing the EEL and the 
SP radiation in three-dimensional photonic crystals com- 
posed of spheresi 21 i 22 ' 23 i 24 

Let's assume that a running charged particle passes 
outside a finite-thick photonic crystal composed of a pe- 
riodic array of the Aluminum cylinders. The particle runs 
with distance s from the boundary normal to the T — X 
direction. The schematic illustration of the system un- 
der study is shown in Figd In this case the induced 
radiation field in the outer region of the photonic crystal 
including the particle trajectory is given by 

E i " d (x,c) = ^ f ^i e *K,; *Q_ +(Mlo)e+5 (39) 



K /s = (k x + h,±ih,k z ), 

— = k x + h , 

v 

lh = yjq% - (k x + h) 2 - k*, 



(40) 
(41) 

(42) 



where Q±±{h,h'), which has spatial tensor index, is 
the scattering matrix of the photonic crystal^ k x is 
the Bloch momentum in the irreducible surface Brillouin 
zone associated with the boundary and h(= 2n7i/a) is a 
reciprocal lattice vector. In the opposite outer region of 
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+e 
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FIG. 7: A charged particle runs with distance s from the 
boundary of the finite-thick photonic crystal, which is com- 
posed of a square array of cylinders with pitch a. The po- 
larization vector of the evanescent wave accompanied by the 
charged particle is denoted by e . 



the photonic crystal the induced radiation field (upper 
SP radiation in Fig[7|) is given by 

E ind (x,^) = E J ^e l <- x Q ++ (Mo) e + . (43) 



Quantitatively, the EEL per unit length of the particle 
trajectory is expressed by the following equation: 

P el (u,k z ) = -^ Q e 2 ioe- 2 ^ s 

x| 7ho |Im[(e+)tQ_ + (/ lo ,/ lo ) e +] . (44) 

On the other hand the SP radiation spectrum per unit 
length is given by 



P em (u),k z ) = \^e 2 ue 2|7fe o 
o 

x £ lh {\Q ++ {hM)e + \ 2 - 

/ifEopcn 



|Q_ + (Mo)e + n,45) 



where the summation is taken over the open diffraction 
channels. 

If the cylinders are loss-less in the photonic crystal, 
again P e i(ui,k z ) — P em (ui,k z ). This can easily be con- 
firmed by considering flux conservation of through the 
photonic crystal: 

]T J h (\Q++{h, h )e+\ 2 + \Q-+(h, h )e+\ 2 ) 

/tfEopcn 

= 2| 7h0 |Im[(e+)tQ_ + (/ l0 ,/ l0 )e+] (46) 

Here, we study two kinds of metallic photonic crys- 
tals. One is a dilute photonic crystal that is composed 
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FIG. 8: The photonic band structure of the square lattice 
of the Aluminum cylinders at k z — was projected on the 
surface Brillouin zone associated with the boundary normal 
to the r — X direction. In (a) the cylinders are well sep- 
arated(lattice constant a — 4r=5[nm]), whereas in (b) the 
cylinders nearly touch(a = 2.16r=2.7[nm]). 



of the square array of the Aluminum cylinders with the 
lattice constant a = 4r = 5[nm]. The other is its dense 
version which corresponds to the cases (c) and (d) of 
the coupled two cylinders treated in the last section 
(a = 2.16r = 2.7[nm]). The photonic band structures 
of these photonic crystals at k z — projected on the sur- 
face Brillouin zone associated with the boundary normal 
to the r — X direction are shown in Fig[H] The band 
structures were calculated by using the two-dimensional 
layer-KKRO method taking l max = 5 and 12 in the di- 
lute and dense photonic crystals, respectively. Here, we 
dropped the band diagram of the S(TM)-polarization, 
because it is not relevant to our problem. However, at 
non-zero k z we must take account of both the S and P 
polarizations. 

In both the photonic crystals the plasma frequency of 
the cylinder is much smaller than the lattice scale, so 
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that the photonic band structure is very close to that of 
the empty lattice at high frequencies. However, below 
u> = ujp many flat bands which characterize the metal- 
lic photonic crystals appear. These bands are generally 
anisotropic, reflecting the C± v symmetry of the square 
lattice, and have a singular state density. In the di- 
lute photonic crystal, many bands are concentrated near 
uj = u> p /\/2. This indicates that they are merely a tight- 
binding coupling of the SPP of the isolated cylinder. In 
principle, we can find infinite numbers of the flat bands 
at uj — uj p /^/2 : which are quite difficult to distinguish. 
On the other hand, in the dense photonic crystal the flat 
bands are diverse in frequency, whereas their center is 
still at uj = u>p/\/2. Some of the bands are originated 
from the cavity mode localized in the groove of the two 
cylinders. However, most flat bands are considered to be 
related with the SPP of an Aluminum cylinder. 

In FiglHl the dispersion lines of the charged particle 
(k x — lo/v + ho) at v = 0.4c and 0.9c as well as the 
light line (uj = ck x ) were also plotted. In the frequency 
range concerned only the line with ho = is relevant, be- 
cause the threshold of the SP radiation occurring along 
the Umklapp shifted line of ho — 2?:/ a is rather high 
(£_ ~ 0.286 and 0.474 in Eq.l^) for v = 0.4c and 0.9c, 
respectively). As a result, the EEL is caused solely by 
the absorption in the frequency range concerned. When 
the dispersion line meets the shaded region of the pro- 
jected band diagram, the charged particle can excite an 
eigenmode in the photonic crystal and thus causes an 
enhanced absorption loss in EEL. Strictly speaking, the 
projected band scheme should be used to understand 
the feature of the photonic crystal with infinite thick- 
ness along the F — X direction. Since we are considering 
a finite-thick photonic crystal, the shaded region in Fig|S] 
must be regarded as a set of the dispersion curves of the 
eigenmodes in the finite-thick photonic crystal. Appar- 
ently, as the thickness increases, the dispersion curves 
fills up with the shaded region. 

FigElshows P e i(co, k z ) with k z — of the two photonic 
crystals, varying number of layers. The velocity of the 
charged particle was taken to 0.4c and the impact param- 
eter s was taken to zero. Concerning the dilute photonic 
crystal, the EEL spectrum has the double peaks near 
uj = uip/y/2 in the mono-layer case. This feature already 
appeared in the result of the two separated cylinders (see 
Fig^J, where an asymmetry of the loss peak is observed. 
As the number of layers increases, the double peaks dis- 
appear and the spectrum converges to a certain function 
which has single peak near uj — uj p /V2. The converged 
spectrum is not so far from the EEL spectrum of the 
isolated cylinder. These features are consistent with the 
numerical results on the projected band structure (Fig[SJ): 
The dispersion line of the particle with v = 0.4c meets 
only the flat bands near uj = ui p /y/2. We can infer that 
the single peak in the EEL spectrum is caused by the 
broadening of the flat bands owing to the non-zero imag- 
inary part in e a . 

As for the dense photonic crystal, there are several 
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FIG. 9: The EEL spectrum of the photonic crystals at k z — 0, 
varying the number of layers. The velocity of the charged 
particle was taken to v = 0.4c. In (a) the cylinders are well 
separated (lattice constant a = 4r), whereas in (b) the cylin- 
ders nearly touch(o = 2.16r). The trajectory of the charged 
particle is just on the boundary of the photonic crystal (i.e. 
s = 0). 



loss peaks whose positions change as the number of lay- 
ers increases. Compared with the case of coupled two 
cylinders, the peak positions of the EEL spectrum in the 
photonic crystal are well correlated with those of case 
(d) of Fig@J. In particular, the two peaks at uj ~ 0.55o;p 
and 0.82w p at N — 32 are of reminiscences of those in 
case (d), and the corresponding flat bands, which have 
relatively large widths in frequency, can be observed in 
FigJSfb). Again, above uja/2wc ~ 0.01 the EEL spectrum 
converges to a certain function with increasing number 
of layers, though the convergence progresses slowly com- 
pared with that in the dilute photonic crystal. A remark- 
able feature in this case appears below uja/2-Kc = 0.01, 
where a frequency shift of a small loss peak is observed 
with increasing N . In contrast to the dilute photonic 
crystal, in such low frequency region the dispersion line 
lo = vk x exists within the shaded region of the lowest 
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FIG. 10: The EEL spectrum of the dilute photonic crystal 
having 256 layers at k z — 0. The velocity of the charged 
particle was taken to 0.4c (dashed line) and 0.9c (solid line). 



band even at v = 0.4c. This band does not originate 
from the SPP modes, and thus the loss peaks found be- 
low ujajluc is different in feature from that by the SPP 
bands. 

The effects of the lowest band in the EEL spectrum 
can be clearly demonstrated in the dilute photonic crys- 
tal with large thickness, using impinging a charged par- 
ticle impinging with such high-speed that its dispersion 
line to = vk x is in the lowest band in FigEl (a). FiglTOI 
shows the two EEL spectra of v = 0.4c and 0.9c in the di- 
lute photonic crystal having 256 layers. As can be seen in 
the figure, the lowest band causes very sharp loss peaks 
whose positions are distributed below aja/2nc ~ 0.025. 
In order to understand this feature, we must remark that 
outside the light cone the band becomes the set of the 
guided modes in the corresponding finite-thick photonic 
crystal. Moreover, as mentioned later, we can reasonably 
introduce an effective dielectric function e c g(w), which is 
very close to that of Maxwell-Garnett^ to the photonic 
crystal under consideration. The effective dielectric func- 
tion of Maxwell-Garnett is given by 



e cff ( LJ ) 



Sb 1 



£6 



2/a 
I- fa 



Sb 



(47) 
(48) 



/ being the filling ration of the cylinders. Using this ef- 
fective dielectric function, the dispersion relation of the 
guided modes in the (loss-less) effective medium is deter- 
mined by 



1 



- 7 / £b + 7 yRc( £ ^) 

7/ eb + 7 '/Re( £o M ff G ) 



(-) a Re(e& G )-** 



exp(2i 7 d) = 0, (49) 
(50) 
(51) 



d being the thickness of the photonic crystal. By im- 
posing the matching condition of frequency u) and wave 
vector k x {= oj/v), the above equation has a sequence of 
solutions, which agree with the positions of the sharp loss 
peaks of v — 0.9c in FigllOlfairlv well. 

The convergence of P e i(uj,k z ) is a direct consequence 

of the convergence of the scattering matrix Q |_ itself. 

As was discussed by Botten et al, the converged value 
of Q |_ gives the reflectance of the semi-infinite pho- 
tonic crystal; 39 This also implies that using the converged 
value of Q |_, we can extract the effective dielectric func- 
tion via Fresnel's formula of the interface S-matrix. That 
is, for the P-polarized incident wave, the scattering ma- 
trix Q |_ of the semi- infinite photonic crystal can be re- 
garded as the interface S-matrix between the background 
medium and the effective medium: 



)-+(ho,ho)]f 



l/s b - V/feg 

7/ £b + 7'/£cff : 



7 



(-) 2 £cff -(-) 2 . 

C V 



(52) 



(53) 



Here, k z = was assumed. The effective dielectric func- 
tion £e ff obtained in this way, along with that of Maxwell- 
Garnett for the dense photonic crystal are shown in 
Fig^J The function is not so far from the effective di- 
electric function of Maxwell-Garnett, though some extra 
features at Loa/lnc ~ 0.019 and 0.026 are observed. In 
the next section we will see that the EEL spectrum in the 
photonic crystal, when the charged particle runs inside 
the photonic crystal, is well reproduced with the effective 
dielectric function having the extra features. As for the 
dilute photonic crystal, our effective dielectric function is 
very close to that of Maxwell-Garnett. 



ELECTRON RUNNING INSIDE PHOTONIC 
CRYSTAL 



When a charged particle runs inside the photonic crys- 
tals, the induced radiation field is rather involved owing 
to the multiple-scattering among the layers above and 
below the trajectory. However, the scattering matrix 
formalism is readily adapted to the case as long as the 
particle does not penetrate any cylinders in the photonic 
crystals. A schematic illustration of the system under 
consideration is shown in Fig^J In this case the induced 
radiation field reacting back to the charged particle is de- 
termined as 



E ind (x,^) =^( a/l e J <- x + b 



he 



iK7 



(54) 



a h = (1 - Q^Ql+y'Q^ie- + Ql + e+), (55) 
b fc = (1 - Q%Qi_r 1 Ql + (e + + Q d + _e-), (56) 

where Q±± is the scattering matrix of the upper(lower) 
layers above(below) the trajectory. Beside, the Fourier 
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FIG. 11: (a) The effective dielectric function for the P- 
polarized light in the dense photonic crystal(a = 2.16r). The 
case k z = was assumed, (b) The effective dielectric function 
of Maxwell-Garnett in the dense photonic crystal. 
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FIG. 12: A charged particle runs inside the photonic crystal 
with an equal distance from the upper and lower nearest layers 
of the trajectory. The plane wave coefficients of the induced 
wave in the void stripe including the trajectory are denoted 
by a and b. 



coefficients of the upper and lower transmitted wave, de- 
noted by tj, is also obtained as 

*ii = E h ')(e +5 Vho + a »'), (57) 

h' 

*fc =E<3--( ft .^)(e"^'fto +bv). (58) 

h' 

Therefore, the EEL and SP radiation spectra per unit 
length becomes 

P ei {uj,k z ) = ^ e 2 u;|7 fco |Im (a£ ■ e" - b ho • (e + )1$9) 
P sp (u,k z ) = ^ e 2 uj -rh(K\ 2 + K\ 2 ). (60) 

/iGopcn 

Again, the flux conservation in a loss-less photonic crystal 
leads 

E 7,(|t^| 2 + |t,-| 2 ) 

/iGopcn 

= -2| 7h0 |Im(a^ -e--b h0 •(€+)*), (61) 

which implies that the EEL is equal to the SP radiation 
emission. 

At low frequencies we may expect that the EEL in a 
photonic crystal is somehow approximated by that of a 
lossy effective homogeneous medium. As is known well, 
the relativistic EEL in such a medium with permittivity 
e cff is given by 



Pei(u,k z ) = i^ e 2 wRe 



I f x ( £ )2 1 



7 



V ScS 



(62) 



per unit length. When e g is real and the condition 
v > c/y/s is satisfied, the above equation is equal to 
the Cerenkov loss. Otherwise, Ea. (|62() can be regarded 
as the EEL by the absorption. In a homogeneous metal 
the bulk plasmon dominantly contributes to the EEL, 
because of the factor l/e e ff- 

FigED shows the EEL spectra of the dense photonic 
crystal and its simulation using the effective dielectric 
function obtained by Eqs 115 2[l and (|53|) . Here, the 
charged particle runs with velocity 0.4c between the 32th 
and 33th layer of the 64-layer thick slab of the dense pho- 
tonic crystal. 

One can observe that the frequency of the main loss 
peak is larger than u> — 0.9uj p , which is close to the bulk 
plasmon frequency. This is quite reasonable, taking ac- 
count that the dense photonic crystal has a large filling 
ratio (~ 67%) and thus is close to the bulk metal of Alu- 
minum. Beside the main loss peak, two small peaks are 
observed in the EEL spectrum. Using the homogeneous 
medium approximation with the effective dielectric func- 
tion obtained in the previous section the above features 
are well reproduced. It should be emphasized that the 
small two peaks in the dense photonic crystal can not be 
explained with the effective medium theory of Maxwell- 
Garnett based on e^ G given by Ea. (|47fl . 
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FIG. 13: The EEL spectrum in the dense photonic crystal 
having 64 layers along the V — X direction is shown. The EEL 
spectrum in the effective medium whose dielectric function is 
given by Fig llll is also plotted. The velocity of the charged 
particle is taken to 0.4c. 

Regarding the EEL spectrum in the dilute photonic 
crystal, our effective dielectric function as well as the 
Maxwell-Garnett approximation reproduce the spectrum 
having a single peak near ui — lu p /^/2 fairly well. 

VI. SUMMARY 

In this paper we have presented a fully-relativistic anal- 
ysis of the EEL and the induced radiation emission in 
various spatial arrangements of metallic cylinders by us- 
ing the multiple scattering method and the layer-KKRO 
method. In an isolated metallic cylinder with a nanoscale 
diameter we showed that the EEL is dominated by the 
absorption rather than the induced radiation emission. 
Thus, the efficiency of converting the kinetic energy of 
the charged particle to the radiation emission is very low. 
In the two identical metallic cylinders a variety of EEL 
peaks appear. Some of them are attributed to the cav- 
ity mode localized in the groove between the cylinders. 
Such a cavity mode as well as the SPP modes become the 
seed of the flat bands in a dense periodic arrangement of 
the metallic cylinders. After presenting a mathematical 
description of the EEL and the SP radiation emission 
in two-dimensional photonic crystals composed of cylin- 



ders, we showed the numerical results of the EEL spectra 
in both dilute and dense periodic arrays of the metallic 
cylinders. In the dilute photonic crystal the EEL spec- 
trum has a simple structure. The spectrum has a single 
peak near cj — uj p /^/2 and is not so far from the EEL 
spectrum in the isolated cylinder. However, when a high 
speed charged particle passes near the photonic crystal, 
a sequence of very sharp loss peaks, which comes from 
the lowest photonic band guided in the finite-thick pho- 
tonic crystal, are observed. The peaks are comparable in 
magnitude with that by the SPP bands. On the other 
hand in the dense photonic crystal the EEL spectrum is 
very complicated reflecting the photonic band structure, 
though a good correspondence to the EEL spectrum in 
the almost touched two cylinders is observed. In both the 
photonic crystals the effective dielectric functions, which 
are obtained with the reflectance of the corresponding 
semi-infinite photonic crystals, fairly reproduce the EEL 
spectra when the charged particle runs inside the pho- 
tonic crystals. 

In this paper we have restricted ourselves to various 
arrays of Aluminum cylinders with a diameter of a few 
nanometers, bearing carbon nano-tube arrays in mind. 
Since the plasma wavelength of Aluminum is much larger 
than the above scale, a metallic photonic crystal com- 
posed of the cylinders behaves as if it has an effective 
dielectric function reflecting the coupled SPPs, in the 
frequency range concerned. As a result, the SP radiation 
from the photonic crystal is completely absent. However, 
it is of great importance to study the EEL and the SP 
radiation emission spectra when the lattice constant is 
comparable with or larger than the plasma wavelength 
of the constituent cylinders. In this case the flat bands 
of coupled SPPs appear in uia/2iTc ~ 0(1), so that the 
EEL is caused both by the absorption and by the SP ra- 
diation. How these bands as well as usual photonic bands 
affects these spectra is the main theme of the paper II. 
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